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Abstract 

The van-der- Waals version of the second virial coefficient is not far from being exact 
if the model parameters are appropriately chosen. It is shown how the van-der- Waals 
resemblance originates from the interplay of thermal averaging and superposition of 
scattering phase shift contributions. The derivation of the two parameters from the 
quantum virial coefficient reveals a fermion-boson symmetry in non-ideal quantum 
gases. Numerical details are worked out for the Helium quantum gases. 

Key words: van der Waals model, Kinetic theory, Quasi-particle methods, 
fermion-boson symmetry 

PACS: 05.20.Dd, 05.30.-d, 05.30.Ch, 31.15.Lc, 51.10.+y 



1 Occupation number statistics and the van-der- Waals model 



Occupation number statistics for N = nfl non-interacting distinguishable 
quantum particles (Boltzmann statistics) in a volume Q yields the entropy 
density functional 

8 = |/dV®v ? (i-h^) (i) 

where vp is the average occupation number of a single-particle energy eigen- 
state. These eigenstates are enumerated by corresponding points p in momen- 
tum space, the density of which is p(p). ks denotes Boltzmann's constant. For 
the more general case of indistinguishable quasi-particles one may consider 
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eq. (22) of [T], where the Fermi-Bose functional is stated. Nevertheless, this 
simplifies to our eqs. (pQ) or <^ if v$ <C 1 which is the case that we discuss in 
this article. 

A hard-core like repulsive interaction is taken into account by the van-der- 
Waals ansatz 

p(p) = (n- Nb)/(2irh) 3 . (2) 

Here the "single-particle volume" b is our first van-der-Waals parameter. We 
consider elementary cells of volume 

- — - ( 27r ^ 3 (3) 
e p{p) 1 — nb 

in six-dimensional phase space (//-space) and the one-particle distribution 
function 

f(P) = VpjVet . (4) 

The entropy density 

s = k B J d 3 pfip) (1 - Hfip)v ee )) (5) 

is then to be maximized as a functional of / subject to the constraints of fixed 
particle density 

n = J d 3 pfip) (6) 

and fixed energy density 

u = Jd 3 pf(p)e f . (7) 

According to the second van-der-Waals ansatz, each particle has its classi- 
cal kinetic energy and is in the potential field of interaction with the other 
particles, i.e. 

ep=- an , (8) 

2m 

where a is the second van-der-Waals parameter. Clearly, the treatment of 
correlations is incomplete in this model. 

The energy density now becomes 

2 

u = J d 3 p^-fip)-an 2 . (9) 
From the principle of maximum entropy s and using the temperature definition 

T = [{ds/du) n }- 1 (10) 
one finds / to be the Maxwellian 

fip) = ni2nmk B Tfl 2 exp . (11) 
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This leads to 

• = ^g--(i^)) < 12 > 

with the thermal wavelength 

A = (2irh)/^2%mk B T (13) 

and to 

u 

2 



3 

u = -nk B T — an 2 . (14) 



Then the pressure formula 



results in the van-der-Waals equation of state: 

^eq = " 7 ~ an ■ ( 16 ) 

1 — no 

Because of the insufficient treatment of two-particle correlations, this formula 
will allow a quantitatively satisfying fit for real systems only if n\ 3 is suf- 
ficiently small. Consequently, the van-der-Waals version of the second virial 
coefficient B(T) is a good approximation if the temperature is not too low: 

B(T)*B vdW (T)=b--?-. (17) 

k B l 

The appropriate choice of the parameters a and b is dealt with in the following 
sections. In particular, it could very well be that both parameters depend on 
the fermionic or bosonic nature of the interacting particles. It will turn out 
that this is not the case. 



2 Heuristics of the van-der-Waals parameters 



The van-der-Waals model can be introduced via corresponding approxima- 
tions to the radial distribution function. To this end consider first the average 
potential energy of iV = nfl mutually interacting classical particles: 

W pot = l<Y, V (\ri-?j\> ■ (18) 

With the radial distribution function defined by 

g(\r'-r"\ >= n' 2 < £)<?(f " '- ri)5(r "' - fj) > (19) 
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one gets 

W pot = N-J d 3 rg(r)V(r) . (20) 
Now g(r) has its density expansion 

g(r) = g (r) + ng x {r) + n 2 g 2 (r) + ... . (21) 

The energy of a single classical particle is therefore - apart from higher-order 
density contributions - given by eq. jSJ) with 

1 Z-,8. 



a=-~ I d 6 rV{r)g {r) . (22) 



In view of the classical limit 

2 M (r) = exp j-^j (23) 

with a Lennard- Jones potential (see below, eq. (135]) ). there will be a cut-off 
radius r* such that the approximation 

, , for r < r* 
9o(r) = { (24) 
1 , for r > r* 

with a constant r* is applicable in a considerable range of temperature. This 
eventually fixes the parameter a as 

a=--f d 3 rV(r) . (25) 

2 J r>r* 

On the other hand, both the parameters a and b may be introduced by first 
expressing the second virial coefficient in terms of go [2] , 

B(T)= l -jtfr(l-g {r)) (26) 
and employing the closer approximation 

, for r < r* 

"""' , = i , w , ■ (27) 

1 - , for r > r* 

Comparison with the van-der-Waals version of -B(T) (eq. (|17p ) then yields a 
as given by eq. (T2"5l) and 

& = y^ 3 . (28) 

The closer approximation of go ~ if inserted in eq. (|22l) - would cause a slight 
dependence of a on temperature. This must be negligible for the van-der-Waals 
model to be acceptable. 
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In the next section, formulae ffTTj) . ff25l) . and (1281) will be substantiated by 
numerical analysis of the exact quantum mechanical virial coefficient, and the 
cut-off radius r* will be determined. 



3 The exact virial coefficient and its van-der-Waals limit 



3. 1 Substantiation of the model 



The exact theory [3] for boson or fermion gases with their two-particle interac- 
tion having, possibly, bound state energies Ei gives the second virial coefficient 
as a sum of four terms: 

B(T) = T 2- 5/2 X 3 - 2 3 / 2 A 3 J2 e- E * /kBT + « G± » + << /V >> (29) 

with the quantities <C G± ^> and <C F± ^> being explained below (eqs. fl37|) . 
(j3HD ) arid with the upper (lower) sign valid for bosons (fermions). The dou- 
ble bracket is our notation for the thermal average of momentum dependent 
functions, e.g. 

POO 

< $ >= / dpw(p) $(p) , (30) 
Jo 

with the thermal weight function 

w(p) = 4n(irmk B T) ^ 2 p 2 ex V { f— 1 . (31) 

[ mk B T J 

Evidently, formula (1TT1) is justified if F± 3> and <C G± ^> prove to be 
practically constant, i.e. -C F± 3>= a and <C G± ^>= 6 in a relevant range of 
temperature where the other contributions are negligible. This is indeed the 
case as will be shown in the following. Moreover, a and b thus defined will 
exhibit a new kind of fermion-boson symmetry in that they are independent 
of the specific quantum statistics. 

The functions F± and G± may be expressed in terms of the properly (anti-) 
symmetrized momentum representation of the two-particle operator 

T(z) = V- V-^—V , T'{z) = ^T{z) , (32) 
H — z dz 

with H = Hkin + V being the Hamiltonian of relative motion: 
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F ± (p) = -h2nh) 3 ^(<p \T ± (E p + ie> \p>) , E p = V - , (33) 
2 m 



G±(p) = -(2nhy J d 3 q5(E p -E g )- (34) 

•9 (< v \T±{E q + ie)\q>< q \T^{E q + ie)\p>) . 

Our graphics Fig. [I] for F± and G± rely on the numerical evaluation for bosons 
( 4 He atoms) interacting via a Lennard- Jones potential lacking bound states 
|4] and fermions (same mass and same interaction as 4 He): 



V{r) = W 



a\ 12 (a 



r 



V /k B = 10.22K, a = 2.56A . (35) 



The (anti-) symmetrized T-matrix is - up to a multiplicative constant - noth- 
ing else but the scattering amplitude 

f±ip-i @) — — ^ 2 mh < p \T±(E P + ie)\q > , \p\ = \q\ , p ■ q = p 2 cos 9 . (36) 

An alternative representation for F± and G± can therefore be given in terms 
of scattering phase shifts 5e which complies with the Beth-Uhlenbeck result 
for B(T) 06]: 

F ± = —f±(p, 0) = —^J2 *W + 1) sin25,(p) (37) 
m m 2p g 

G±(p) = -Avh%Y. ± W+ l)sin 2 [5,(p)]^|M . (38) 
p z £ op 

The summation runs over even I for bosons and odd i for fermions. 

A remark on the units and dependencies of F±, G± will be fitting here. Let 
p = hk and choose f = r/a as the dimensionless radial coordinate. Then the 
dimensionless version of the radial wave equation with eigenvalue E p becomes 

u't(r) + {Re (r 12 - f ~ 6 ) + £(£ + l)f~ 2 } u(r) = (ka) 2 u £ (f) (39) 

where 

Re = (40) 

is the "Reynolds number". Given the appropriate behavior of ui for f — > 0, 
the asymptotic behavior for f — > oo exhibits the phase shift Sf 

£ 

Uiif) — > sin(kar — -tt + 5f) . (41) 

Thus the only dependency of 5e is on Re and ka. Consequently, with some 
functions <p± = y?±(Re, ka) and j± = 7±(Re, ka) 

„ Anh 2 a , . n (j9+(Re, ka) , . 

F± = <^± Re, fox = 16nV a 3 , 42 

m Re 
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kd kd 



Fig. 1. Functions F± and G± for bosons ( 4 He atoms) interacting via the 
Lennard-Jones potential (|35p and fermions with the same mass and same inter- 
action. The thermal weight function is given for two temperatures. The horizontal 
line on the l.h.s. marks the value of a as given by (|59p and on the r.h.s. the value 
of b as given by ([551) . 

G± = 47rt 3 7 ± (Re,fca) . (43) 

In each of our graphics the boson and the fermion function refer to the same 
value of Re. With the potential data of eq. ([351) and with m the mass of 4 He 
we obtain Re m 22.1. 



The weight function w = w(p) (eq. ([3TI) ) assumes its maximum at p Tl 
hk max with 



/ lOORe \ 1/2 / T x 1 / 2 
k — 0= \A^22) \i00k) (44) 

and it is easily seen that above T = 100 K the weight functions samples rather 
large values of ka iko ^ 10). For high values of ka Fig. [1] shows that both 
F + and F_ oscillate about a common approximately constant value. This 
reflects the nearly hard-core likeness of the repulsive part of the Lennard- 
Jones potential. Let us therefore consider, for comparison, a pure hard-core 
repulsion with radius a. 

In this special case 

tan5 £ = j e (ka)/y e (ka) (45) 
holds, with jg and yt denoting the spherical Bessel functions 

sinz / d \ n cos 2; fi . 

, y e (z) = z n [ — ) . 46 

z \dz J z 

Invoking 

tan 5p 

sin 25 1 = 2- 47 

l + tan 2 5 £ v y 
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and the asymptotic behavior 



hm(j^)+^)) = l/z 2 



(4£ 



one arrives at 



lasy 



A-Kt) 2 cr 2 k 00 

E(-!) m (^ + 1) * ^a) ye(ka) 



(49) 



Now - as a marginal case of formula 10.1.46 in [7] 



£(-l) £+1 (2£ + l)j i (z)y l (z) 



cos2z 
2a 



(50) 



and consequently 



he 



?hc\ 



2nh 2 a 



lasy 



m 



cos(2A;cr) 



(51) 



This is the asymptotic (ka — > 00) result for the hard-core system. It is com- 
pared with (F_ — F + ) for the Lennard- Jones system in Fig. [2J The behavior 
is very similar both in terms of amplitude and frequency. A slight decrease of 
a with increasing k in formula fl5Tl) would still improve the agreement. This 
reflects the fact that the Lennard- Jones potential appears the softer the higher 
the particles' energy is. 
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Fig. 2. Comparison of (F- — F+) for the Lennard-Jones system (solid curve) with 
the asymptotic result for the hard core system (dashed curve). 

As for G±, an oscillation about a common constant value is once again seen, 
compare the r.h.s. of Fig. [IJ In contrast to the case of F±, however, the oscil- 
latory amplitude of the difference is clearly decreasing. Not surprisingly, this 
feature can again be derived analytically for a pure hard-core repulsion with 
radius a. 
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To this end, G± (eq. (!38l) ) is first rewritten as 



(52) 



and then employed for the hard-core system. With the aid of eqs. (BSD, (HT|) . 
and (148]) the asymptotic (/ca — > oo) tail of (G+ c — G^. c ) is found to be 



(G h + C - G hc " 



lasy 



27ra 3 ^(-l)^(2£+l)j,(Hy;(^) • 



(53) 



After replacing (for ka — > oo) y' s by je, one can apply formula 10.1.51 in [7], 

sin 2^; 



£ (-1)^ + 1^ 

£=0 



and hence 



2^ 

sin 2/cer 



(Cr + - G_ ) ass/ = 27TCT 



(54) 



(55) 



The difference (G+ — G_) for the Lennard- Jones system is compared with this 
result in Fig. [3j Again, the behavior is very similar both in terms of amplitude 
and frequency, also the phase difference increases only slightly. 



"Lennard-Jones 
hard core, asymptotic 




Fig. 3. Comparison of (G_ — G+) for the Lennard- Jones system (solid curve) with 
the asymptotic result for the hard core system (dashed curve). 



3.2 The cut-off radius 



Instead of eq. ff271) a continuous ansatz for go(r) may be used: 

, , for r <* 

9o(r) = { , y(r)-V(r.)wp(-«*=?*) (56) 
-a^f) + ^ *- , for r > r* 
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with the potential V(r) according to eq. (|35|) . Equation (|2"T|) is reproduced for 
a — > oo . The values of r*/<r and a follow from the van-der-Waals condition 

'-if^y^ 1 -^ < 57 > 

where 

b = 0.4trj 3 (58) 

and 

a = = -5- • V a s 59 

m Re 

is estimated in view of Fig. [1] (horizontal lines). Then 



where x is the solution of 



r* = x 1/3 a 



x = ^- + 0.3a; 2 + -x 3 - 0.3a; 4 (61) 
Re 3 

with 0.2 x > 1/3. The parameter a follows from 

0.2a; _ 1/3 = - + — + • (62) 
a cr cr 

Then a has the following representation in terms of the Reynolds number Re 
and the Lennard- Jones parameters V and a: 

47r"iV 3 ( 1 3 



(60) 



a = [x - -x 3 ) , (63) 

where x = a;(Re) is that one among the solutions of eq. flBTj) . which allows 
a positive a in eq. (I62lh For Re= 22.1 we obtain x = 1.88 and a = 25.6. 
Therefore, in our ansatz fl56l) go( r ) is close to its limit for a — > 00, which is 
given in eq. (T2~7|) . 



Summary 



The summary of our elaboration is that for a large region of temperature 
the functions F± ^>= a and G± ^>= 6 may be nearly considered 
as constants, i.e. not depending on temperature, see Fig. [TJ Moreover, the 
respective kind of statistics (Bose-Einstein or Fermi-Dirac) does not matter. 
This constitutes a fermion-boson symmetry in non-ideal quantum gases. For 
ideal quantum gases such symmetries are known for about ten years [8119] . 

The slight residual dependence of a and b on temperature (cf. Fig. [[]) reflects 
the fact that these are parameters not of an exact but of a model theory. 
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